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The aim of this paper is to extend linear quantum dynamical network theory to include static 
Bogoliubov components (such as squeezers). Within this integrated quantum network theory we 
provide general methods for cascade or series connections, as well as feedback interconnections using 
linear fractional transformations. In addition, we define input-output maps and transfer functions 
for representing components and describing convergence. We also discuss the underlying group 
structure in this theory arising from series interconnection. Several examples illustrate the theory. 
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I. INTRODUCTION 

The aim of this paper is to develop systematic meth- 
ods for describing and manipulating a class of quantum 
feedback networks. Quantum networks are important 
because of the fundamental role they play in quantum 
technology, and quantum information, computing, and 
control in particular. Knowledge of how to interconnect 
quantum (and classical) components in a network is a 
pre-requisite for feedback control. 

A quantum feedback network (QFN) Q, Q is a phys- 
ical network whose nodes (components) are open quan- 
tum systems and whose branches (wires) are quantum 
fields. Quantum dynamical components have linear rela- 
tions between their input and output fields, but in general 
their internal physical variables may evolve nonlinearly. 
Furthermore, delays may be present due to the non-zero 
physical length of the branches and the finite speed of 
light. The network theory developed in Q, Q is based 
on quantum stochastic models of open quantum systems 
(Hudson-Parthasarathy [l^l, Gardiner-CoUett [1]), and 
provides methods and tools for QFN modeling including 
series (or cascade) connections, and feedback loops. The 
series connection defines a group operation in the class of 
open quantum dynamical systems. In this paper we will 
consider the subclass of networks consisting of dynamical 
components whose internal variables evolve linearly. In 
[Hi the theory was presented for components based on 
unitary transformations such as beamsplitters and phase 
shift modulators, but does not include static components 
that require an external source of quanta for their oper- 
ation as such amplifiers and squeezers. 

In quantum optics one encounters the class of quantum 
linear networks with components implementing static lin- 
ear transformations, called Bogoliubov transformations, 
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see ■ Examples of Bogoliubov components include de- 
vices capable of creating squeezed states of field ( "squeez- 
ers") out of a vacuum, a non-unitary process. Series 
connection of static Bogoluibov components is given by 
matrix multiplication of representations of the Bogoli- 
ubov transformations, and defines a group product for 
this class of systems. 

These two classes of quantum networks, which are 
characterized by the nature of their components, are dis- 
tinct, but have elements in common. The beamsplitter 
and phase shift are components in both classes. However, 
the squeezer does not belong to the class of open dy- 
namical quantum^ systems (in the framework of Hudson- 
Parthasarathy [lOj] and Gardiner-CoUett Q), though it 
may be approximated by systems that are in this class. 

The purpose of this paper is merge together these two 
classes of linear networks in a unified, multivariable al- 
gebraic framework. By 'multivariable' we mean that the 
framework allows for systems comprised of multiple os- 
cillator modes and multiple field channels; accordingly, a 
vector-matrix notation is used. The new class of linear 
QFNs (LQFNs) we consider in this paper are therefore 
assembled from components of the following two types: 
(i) dynamical components, with linear evolution of phys- 
ical variables, and (ii) static components characterized 
by Bogoliubov transformations. An example of such a 
network is shown in Figure [H [l^ , Q . 

To this end, we consider linear open quantum compo- 
nents that are in general a series connection of a linear 
dynamical part, and a static Bogoliubov part. We de- 
fine series connections of these components, and extend 
linear fractional transformation (LFT) methods for de- 
scribing feedback loops. The series connection defines a 
group structure for this new class of systems, which in- 
cludes the linear dynamic and static Bogoliubov classes 
as subgroups. This group structure is interesting from a 
physical as well as a systems and control theoretic point 
of view. 

Open quantum systems have a natural input-output 
structure. We define and make use of input-output maps 
for our class of linear open quantum systems, and discuss 



2 

















sq 












yi 






beamsplitter 






(«) 


U2 


delay r 








V2 



cavity 



FIG. 1: Linear quantum feedback network (LQFN) consist- 
ing of a beamsplitter, squeezer, and cavity. The time delay 
around the optical loop is r. 



convergence of systems in these terms. This input-output 
notion of convergence is important for apphcations, and 
is weaker than stronger notions of convergence involving 
all system variables. 

We begin in section |TT] by describing Bogoliubov trans- 
formations, which is followed in section IIIII by a dis- 
cussion of open linear dynamical models of Hudson- 
Parthasarathy type. In section |V] we discuss quantum 
components involving Bogoliubov transformations, both 
static and dynamic. Linear quantum feedback networks 
are described in section IVII Several examples are dis- 
cussed in sections |IV] and IVIII 
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where we introduce the 2r x 2k doubled-up matrix 
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We note that A{E_,E+y = A{El,El), and when 
the dimensions are compatible, A{E-, E+)A{F-, F-i-) = 
A{E^F_ +E+F*, E_F+ + E+F*). In the examples we 
consider, the linear transformations will be between vec- 
tors of equal dimensions and so the matrices E±, etc., 
will be square. 

For a 2n X 2m matrix X, we define an involution b by 
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where 
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with /„ the n X n identity matrix. When understood, we 
shall often drop the dimension index and just write / and 
J. For doubled-up matrices, we then have 



A{E^,E+)^ = A{Ei,-El) 



(5) 



II. BOGOLIUBOV TRANSFORMATIONS 

In this section we present models for the quantum com- 
ponents considered in this paper. Before this can be 
done, some notation is needed. 

A. Notation 

Let X = (Xjk), j,k = 1, . . . , n, denote a matrix whose 
entries Xjk are operators on a Hilbert space Sj, or are 
complex numbers. We define the matrices 

x* = ix;,), x^^iXk,), x^ = {xi^). 

Here, the asterisk * indicates Hilbert space adjoint or 
complex conjugation. 

For a column vector x of operators of length fc, we shall 
introduced the doubled-up column vector 

of length 2/s, so that ft = [x\x^). 

Given a linear transformation of the form 

y ~ E^x + E+x* 

where x and y are vectors of operators of lengths k and r 
respectively, and E± G ([yxk ^ define the transforma- 
tion y"^ = E'Ulx'^ + E'^x, and in doubled-up notation we 



B. Canonical Commutation Relations 

We consider a collection of m harmonic oscillators, 
whose behavior is characterized by independent annihi- 
lation aj and creation a* operators (j = 1, • • • , m) satis- 
fying the canonical commutation relations [oj, aj^] = 5jk, 
with [aj,afc] = = [a*,aj^]. The commutation relations 
may be written compactly as 

where {Jjk) = Jm is the matrix defined in ([4]). 

Systems consisting of m oscillator modes are equiva- 
lent, for fixed m, and it is convenient to consider just 
the category iS(m) of such systems with representative 
described by column vector a = (ai, . . . , am)^ ■ 

C. The Bogoliubov Matrix Lie Group, Sp(C'") 

Definition A 2m x 2m complex matrix S is said to be 
b-unitary if it is invertible and 

S ~ SS^ ~ l2m- 

The group of Bogoliubov matrices Sp(C™) is the sub- 
group of \? -unitary matrices that are of doubled-up form, 
that is ^ A(S'_ , 5"+) for suitable S- , S+£ C™^" . This 
is also known as the symplectic group flSl L 
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The transformation a' = S-a + S+a^ is called a Bo- 
goliubov transformation for a g S{m). In doubled-up 
notation this takes the simpler form 



Sa. 



(6) 



Note that a' G S(rn) and in particular, the transforma- 
tion preserves the canonical commutation relations. 

A Bogoliubov matrix S S Sp(C™) admits a Shale de- 
composition [13] 

S = A(S'^^t,0)A(coshi?,sinhi?)A(S'i„,0), (7) 

where Sin, S'out arc mx m unitary matrices, and R a real 
diagonal mx m matrix. Note that A(coshi?, sinhi?) = 
expA(0,i?). The middle term in ([Tj) corresponds to 
squeezing, an important characteristic widely exploited in 
applications of quantum optics. To see what this means, 
suppose S = A(coshi?, sinhi?). Define the quadratures 
= ^(a + o,'^) and = ^^^(a — a"^), and similarly for 
a'. Then 



{ay 



{ar 



-""ay. 



(8) 



which shows that if the y quadrature is scaled by less 
than unity, the x quadrature must correspondingly be 
expanded by an amount greater than unity. Also, note 
that the unitary group U{m) of unitary m x m matrices 
can be viewed as a subgroup of Sp(C™) via the corre- 
spondence A{S,0) e Sp(C'") whenever S e U{m). 

The Bogoliubov transformation ([6]) defined by a fixed 
S e Sp(C'") corresponds to the action of a physical de- 
vice acting on a vector a e S{m). By Shale's theorem 
[l^ the Bogoliubov transformation may be unitarily im- 
plemented, that is, there exists a unitary operator U such 
that 

Sa = U*aU. (9) 

D. The Bogoliubov Lie Algebra, sp(C'") 

We remark that the Lie algebra sp(C'^) consists of ma- 
trices —ifl G c2mx2m ^^i^^ g^j.^ doublcd-up form (in 



order to generate doubled-up matrices S 



-in 



) and 



satisfy = f2. The second condition can be written as 
Jf2 — fjt J = 0. We therefore deduce that the infinitesimal 
generators take the form 



-in^ -A(iO_,zO+), 



(10) 



with complex matrices 51_ and having the symmetries 
= il_ and fij = We remark that we may con- 
struct a hcrmitcan operator H on the oscillator Hilbert 
space as 



H = 



1 



1 



Q,/3 = l 



(a>a/3a/3 + i^<apUJaf3 + ^ao^apuj^*), (11) 



where the coefficients are the entries of the matrices 
^± = ('^Q/j) € C™^"'. From the familiar quantum me- 
chanical point of view, H is the Hamiltonian generating 
the canonical transformation ([5]), that is, in ^ we have 
U = e-'". 

It is instructive to look at the m ~ 1 case. Here the 



Hamiltonian is H = uj_a*a -\- ^u^a*'^ 



^Lu'^a with 

real and cij+ = jw+le*^ complex. The corresponding el- 
ement — iil e sp(C) can be written in terms of Pauli 
matrices as as 



-in 



— A(za;_, iujj^) 



where lua 



IUJ4 



or 



UJ+X + itO+y 



The eigenvalues of —iSl are ±-\/^ where 



(12) 



We note that Heisenberg dynamical equations will be 
trigonometric type for ^ < 0, and hyperbolic type for 
C > 0. Let us try and diagonalize the Hamiltonian 
by introducing the Bogoliubov transformation e'^a = 
cosh r a! — sinh r a'* , (this is the original purpose of Bo- 
goliubov transformations!). For ^ < we may choose 

tanh2r = to get H = uj^,Jl 

While for C < we may choose tanh2r = lo 



get H 



1 



a a . 
uj+\ to 
That is, for C < 



wc may diagnolize f2 as S VtS ~ A(tj(,,0) using a Bo- 
goliubov matrix S*, but for C > the best we can do is 
to put into the form A(0,a;![_). 

We say that H is passive if f2 has only real eigenvalues. 
In this case we may find an S in Sp(C) such that S^VlS = 
A(f2'_,0) for some The term "passive" means that 
such Hamiltonians do not describe energy flow into the 
system from an external pumping source, and that the 
dynamical equations are always trigonometric type. 

(For TO = 1, if is passive if and only if the parameter 
C < in as the eigenvalues are i-^^. ) 

The group Sp(C'") is a non-compact group and, in 
fact, is not covered by the exponential mapping on its 
Lie algebra sp(C™). We see this in the case to = 1, 

where the form of e^'^ depends on the sign of C- We 

have e~*^ given respectively by (C > 0) 



* sinh 



BinhVC 

+ Vi 



and (C < 0) 



cos ^/--C, — iuj- 



sin 
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(note that fi^ = in this 

e > -2 
e Sp(C) which do 



we must have tr " 



161 for 



Also e-'" = 1 - in if C 
case). As observed in [l(__ 
so that there exist matrices in S 
not possess a logarithm in 5p(C), for example. S — 
— A(cosh?i, sinhu) = exp(— A(i7r, 0)) cxp(— A(0, — u)). 
In particular, such Bogoliubov transformations are not 
generated by a single Hamiltonian H. The best that can 
be done is to write the unitary t/ in (l9|) as t/ = U\ ■ ■ ■ Uk 
where each Ui has a logarithm in sp(C"'), see 
higher-order cases. 



E. Gaussian States 

A state on S (m) is said to be Gaussian if wc have 

(expi = cxp — — ?l^Fit + iu^aj- , 

where F > 0. For simplicity we consider mean zero states 
(a = 0). In particular, we have that F 
specific form 



(aa^) takes the 



F 



A/t N 



with 



Njk = {a*ak), Mjk = (ojafc) 



(13) 



(14) 



and wc note that N = iV^ and M = . In particular, 
positivity of F implies that > 0. The vacuum state is 
the special state determined by the choice N = {],M = 0, 
for 



/ 




(15) 



For fixed A^ > the choice of M is constrained by the 
requirement that F be positive. For the n = 1 case, A^ 
and M are scalars and the positivity condition is eas- 
ily seen to be A^ > with |Afp < N{N + 1). More 
generally we should have a diagonalization V^NV = 
diag (A^i, • • • , A„) for unitary V in which case we could 
consider new fields a' = Va. Here Nj can be interpreted 
as the average number of quanta in the mode a^- . In gen- 
eral we cannot expect to simultaneously diagonalize A^ 
and M. 



1. Generalized Araki-Woods Representation 

Given a Gaussian state determined by F in equation 
(|13p . we now show that we can construct modes having 
that state through canonical transformations of vacuum 
modes. That is, given a state for which a e S{m) has 
covariancc F given by (|13[) , there exists a 2m x Am matrix 
So such that 



where ao 
and 



ai 

0,2 



G S{m + to) has vacuum statistics. 



5'o'S'o = I- 



(17) 



Indeed, we will construct So — A(£'° , E'l_) for some to x 
2to a matrices E'^, E^. This generalizes a construction 
originally due to Araki and Woods [2l| for non-squeezed 
thermal states, see 123. 1231. 



2. Construction of Araki- Woods Vacuum Representation 

Step 1: Diagonalize N . We may find a unitary matrix 
V e C"^'" such that V'fNV = diag (A^i, • • • , A^„). The 
eigenvalues are assumed to be ordered such that A'^i > 
• • • > > 0. Therefore we can restrict our attention 
to the case where A^ is diagonalized in this way. 

Step 2: Ignore zero eigenvalues. Take the first to-|_ 
eigen values to be strictly positive, with the remaining 
Too ~ m — m_|- to be zero. We respect to the eigen de- 
composition C" = C'"+ ® C™" we decompose F as 



F 



However, we observe that if a positive matrix has a zero 
on a diagonal, then every entry on the corresponding row 
and column must vanish [26j so that actually 



I + N++ 





M++ 


M+o 





/ 


Mo+ 


Moo 






N++ 





A/oV 


Mot, 









F: 



I 





- 


A/++ 


0- 





/ 













N++ 














0. 



Therefore, we can restrict our attention to the case where 
A^ is diagonal and strictly positive, and in particular in- 
vertible. 

Step 3: Explicit Construction. We begin by noting the 
constraint I + N > Mj^M^ which follows from noting 
the positivity of 



" / 




F 


I 


-Af — ' 


t 



















I + N - 



Mj^M'^ 







We then introduce the following matrices: 
X = 

V A 

/aT, • • • , ^/nZ) 



Note that Y = Y^ and from Z = MY^^ we have that 
YZ^ = = M = ZY = ZY^ . These matrices satisfy 
the conditions 




S'o&o 



(16) 



XX^ - YY^ + ZZ^ = I and Y Z'^ = ZY^ 



(18) 
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Now take &i and 62 to be independent (commuting) 
modes in iS (m). We fix the state to be the joint vacuum 
state for both of these modes. Then we may represent a 
as 

a^Xbi+ Yh\ + Zh2\ (19) 
indeed, it is straightforward to check that 

(aa^) = ZY^ = ZY = M. 

Therefore we have obtained the representation (fT6|) 
with Sq = A{E^,El) and 

E°_ = {X 0) , E\ = {Z Y) . 

Property (H?]) follows from (fTS]) . 

III. QUANTUM OPEN LINEAR DYNAMICS 

In this section we consider the general class of open 
linear dynamical models arising from a unitary model 
for the joint system and field. The system will be a col- 
lection S{m) of m harmonic modes with representative 
a = (ai, . . . , amV . 

A. Boson Fields: Vacuum States 



and in the vacuum representation their future-pointing 
(Ito) increments satisfy the quantum Ito table 



X 


dBl dBk 




Sjkdt 


dB] 






We may write the above more compactly as 
(expi / v){t)b{t)dt\ = 

\ Jo I vac 

1 

exp-- / il'' (t) Fys.cu{t) dt. 
^ Jo 

The vacuum state is then the Gaussian state for which 
( b (t) 6t {t') ) = Fvac'5(t - t') and the Ito table may be 

\ / vac 

summarized by 

dBdB'' = F^^cdt. 

Here B = {Bi, B2, ■ ■ ■ , Bn)^ ■ In the vacuum case we may 
also define the counting process 

^Mt)= f b*{r)bu{r)dr, 
Jo 

which may be included in the Ito table (TOj. The addi- 
tional non-trivial products of differentials are 

dAjkdBj = SkidBj,dBjdAki = SjkdBi,dAjkdAii = SkidAj^. 



The open quantum systems to be described below are 
driven by n quantum noise fields (input processes) rep- 
resented by annihilation bj{t) and creation b*{t) opera- 
tors (j = 1, • • • , n) satisfying canonical commutation re- 
lations [b.i{t),blit')] = 6jk5{t~t'), with [bj{t),bk{t')] 0. 
This may be written compactly as 



(20) 



We shall denote the class of m independent input pro- 
cesses by T{n) with representative described by column 
vector b = (&i, . . . , 6„)^. 

The vacuum state for the field is characterized by 



exp I / {u{t)b^t)+u^ {t)b{t))dt 
Jo 



exp ^ 



1 



v) it) u {t) dt. 



It is convenient to introduce the integrated fields 



Jo 



B. Boson Fields: Gaussian Field States 



We may generalize the situation in IIII AI to the case 
where the input fields are in Gaussian states with zero 
mean, but with the correlation functions 

{b*{t)bk{t'))^N,jk S{t-t'), 
{b,{t)bk{t')) ^ M,k S{t~t'), 



(21) 



with N and M as in That is, 

J6(t)6t {t'))^F5{t~t')., 



where F has the same form encountered in the case of a 
finite number of modes in equation p3p . The extended 
Ito table is then 



dBj 
dB\ 



dB,, 



i6jk+Nkj)dt Mjkdt 
Njudt. 



M*^dt 



Generalized Araki-Woods representations for arbi- 
trary Gaussian field states may be obtained based on 
a straightforward lifting of the constructions in III El 

MM. 
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C. Quantum Linear Dynamical Models 



and 



The dynamical behavior of a system comprised of m os- 
cillators, interacting with n input fields (vacuum state) is 
given in terms of the Hudson-Parthasarathy Schrodinger 
equation (in Ito form) 

dU{t) = I {S,j - 5,,) dK,^{t) + ^ dB*{t) 

n 1 1 

- L*S,kdBu{t)-{-YL*L,+iH)dt\u{t),{22) 



for a unitary operator U{t), with U{0) = /. To obtain 
a unitary evolution leading to linear dynamics we must 
take 5* £ to be unitary, H to be of the form encoun- 

tered in while the coupling of the system modes to 
the fields is to be of the form 



(23) 



where C± = {CfJ e C"^™. 

The oscillator variables evolve unitarily aj{t) = 
U*{t)ajU{t), and likewise the output field is i?out(0 = 
U* {t)B{t)U{t). The dynamical equations are 

a{t) = CjS* b*{t) - ClSb{t) + A^a{t) + A+a*{t), 
boutit) = Sbit) + C-a{t)+C+a*{t), (24) 

where 



(25) 



Note that -^(A_-A^) = and -j-{A+ + Al) = n+, 
but in general A_ 7^ A[_ and A^ 7^ A^. Here and below 
differential equations are expressed in terms of quantum 
noise fields, and may be interpreted in the Stratonovich 
or Ito sense, moreover the evolution preserves the com- 
mutation relations of the oscillator variables. 

The linear dynamical equations can be written in 
doubled-up form as 

hit) = A(A_,^+)a(t)-A(C_,C+)^A(5,0)6(i), 
kntit) = A(C_,C+)a(i) + A(5, 0)6(0. (26) 

Let us introduce the doubled-up matrices A = 
A(A_,A+), C = A(C_,C+), and -in = -A{in_,in+), 
then we have the identities 

2Rei,{A) = A + A^ = -C^C, = (l. (27) 

These are readily established by noting 

A(^_ ,A+) + A{A^,A+f ^ A(^_ + Al ,A+-Al) 
= -A{ClC- - CjC*, ClC+ - ClC*) 
= -A(C_,C+)^A(C_,C+), 



A{in_,in+f = A{-in^_,-inl) = -A{in_,in^ 

The dynamical equations can then be recast as 

h{t) = Aa(t) + Bb{t), 
kutit) = Ca{t)+Db{t), 



where 



and 



D = A{S,0), 



B = -C^D, 



A = -^C^C 



in. 



(28) 



(29) 



(30) 



We denote this class of linear Hudson-Parthasarathy 
systems with n input fields and m oscillators by 
C^^ {n,m), and write C^^{n) = Um,C^^ {n,m). Systems 
G G £^^(71,771) may be parameterized in several ways. 
In terms of the scattering matrix, 5", vector of coupling 
operators L and Hamiltonian H we may write 



G= iS,L,H). 



(31) 



Since these physical parameters are determined by the 
matrices given above, wc may also write 



G^is,c,n). 



(32) 



Alternatively, we may use the matrices appearing in 
equations 



G 



'A 


B ' 


c 


D 



(33) 



a notation commonly used in linear systems and control 
theory. We remark that an arbitrary quadruple of ma- 
trices A, B, G, D need not necessarily correspond to a 
quantum open system, [ll| , [T^ . 



D. Stability 

In linear systems theory, a system G of the form ([33]) 
is said to be Hurwitz stable if the matrix A has all eigen- 
values having strictly negative real part. Now if the joint 
system-field state is one in which the inputs are mean- 
zero, then ^(a(t)) = A{a{t)) and so Hurwitz stability 
implies that (a (t)) — > as t ^ cxd. 

When A+ = 0, we have A = A(A_,0) with A^ = 
-^GiG- - in^ and n*_ = n^. Since X^X is non- 
negative definite it is easy to determine whether A- is 
Hurwitz. For instance, it is sufficient to have C_ invert- 
ible. However, expressions such as X^X are indefinite 
due to the presence of the matrix J. There may be non- 
passive contributions to A from both C+ and 

As an illustration let us consider how the eigenvalues 
of A depend on the physical parameters in the simplest 
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case n = 1 = TO. We have seen that the most general 
parameterization is 



C = A(V7^e''^-,V7Te"^+) and A , w+ ) , 
with 7J-, (f)± and ti>_ real and lu^ £ C. In this case 



A 



-A(- (7- 



7+) + iuj_,iuj+) 



which has eigenvalues ^ (7- — 7+) ± where we recall 
the parameter ( — \u!+\'^ — uPl from (fT2|) . The plant is 
Hurwitz if 



1. < and 7_ > 7+; 

2. or C > and ^/C < 5 (7- 



7+). 



In situation 1 the system has a passive Hamiltonian 
and the damping rate is greater that the pumping rate. 
However situation 2 shows that if the damping is suf- 
ficiently large then the system may still be stable even 
if the Hamiltonian is not passive. In general, as one ex- 
pects, stability will depend on the relative flows of energy 
into and out of the system. 



E. Series Connections 

Open linear dynamical systems G\ = {Si, Li, Hi) and 
G2 = (5*2, ^2, H2) in (recall the parameterization 

([3T|l ) may be connected in series by passing the output 
of system Gi into the input of system G2, @, 0, @- 
The system formed from this connection in the zero-delay 
limit is an open system G = G2 < Gi, which in terms of 
the parameters (j31[) is given by 



G2<jGi = {8281, L2 + S2L1, Hi + H2 + lm{LlS2Li}). 

(34) 

We refer to <i as the series product of Gi and G2 and 
gives the cascaded system. The set (n) forms a group 
with respect to the series product, with inverse G~^ ~ 
(5t, -S-tL, -H) (where G = (5, L, H)). 

Given an open system G = (5, G, fl) (now we use the 
parameterization (j3ip in anticipation of later use), it fol- 
lows from properties of the series product that 



G= (/,G,O)<](5,0,0). 



(35) 



This factorization says that an open system with scatter- 
ing 8 is equivalent (in the zero-delay limit) to a dynamic 
open system without scattering (/, G, fl) connected in se- 
ries with a non-dynamic or static open system (S*, 0, 0). 



F. Input-Output Maps 

In classical systems and control theory the input- 
output map is a basic tool, which in the case of linear 



systems may be expressed explicitly in the time and fre- 
quency domains. Input-output maps for the open quan- 
tum linear system may be defined in the same way; in 
terms of the doubled- up parameters p3p . we have 



fcout(0 = Ce^a(0) + SG[t;&in] 



(36) 



where 



Ce^^^-'-^C^Db{r)dr + Db{t). (37) 



Here the input 6i„ is understood as the input field b 
and we shall often use the subscript for emphasis. The 
impulse response associated with the term is 
acit) = -Ce^*&b -f b6{t), from which we have the 
transfer function (the Laplace transform of in 
which s is a complex variable): 



Sg(s) 



" A 






D 



(s) = -C{sl - A)-^&b + b. 

(38) 



Let us introduce the transformed 



e-^*6in (t) dt, 



*6i„ (t) dt and b* [s] = 6i, 



(39) 



that is, &in [s] — 

!^e-^%*{t)dt. ^ 

Adopting a similar convention for the outputs, we then 
obtain an input-output relation of the form &out [s] ~ 
Eg,- (s) 6in [s] + Sg,+ (s) &*n [s], or 



where 



Sg(s) 



&out [s] = ^g(s)V [s] 



Eg,- (s) Eg,+ (s) 
EG,+ {s*f EG,-is*f 



(40) 



(41) 



and we have ignored the initial value contribution of the 
system modes. 

Note that while the transfer function Eg{s) is uniquely 
determined by G, the transfer function does not uniquely 
determine the system G - many systems may have the 
same transfer function. 



IV. EXAMPLES 

A. Annihilation Systems 

A system G = {S,C,n) with G+ = 0, f7+ = 
has dynamics and output relations that depend only 
on the annihilation operators and annihilation fields, 
0, [13 • For a physically motivated reason, since 
neither the Hamiltonian (passive) nor the coupling oper- 
ator of the system contain terms that would require an 
external source of quanta (i.e, a classical pump beam) to 
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implement (this follows from the synthesis theory of [l8| ; 
see also [23, Section 7] for a discussion), they are also 
referred to as passive systems [131 ■ This type of system 
often arises in applications, and includes optical cavities. 
Transfer functions for this class of systems takes a simpler 
form, as we now describe. 

We have A_ = -iclC_ + i^l- and A+ = 0. Then 
the matrices C = A(C_, 0) and A = A(A_, 0) are block 
diagonal and the transfer function takes the form 



Sg,- (s) 
Sg,- {s*f 



(42) 



with 



Sg.-(.s) 





-els' 




s 



(s) = -C-{sI-A.)-^CiS+S. 

(43) 

In this situation, wc have the input-output rela- 
tion 6out(t) = Ce'^*a(0) + ^G[t;b] with ^G[t;b] = 
- jlCe^'^*-''^&Db{r)dr + Db{t). In comparison with 
(f351) and ([ST]) , the output field depends affinely on 6, but 
not the conjugate li^ . 



u+ = |e, e > 0, C- = -/k, and C+ = 0. The Hamilto- 
nian will not be passive, however, the system is stable in 
the sense of Hurwitz if take e < k as we shall do from 
now on. Using ([55]) we find that the doubled-up transfer 
function, in agreement with Q, is 



^dpa(s) 



1 



P{s) 



where P{s) = [s + \tif' — ■ The poles of the transfer 
function therefore occur at the zeros of P, namely s = 
±1 — ^. In the frequency domain, the output field is 

1 . o k2 -He\, _ 1 



6out(s) 



P{s) 



-)b{s) 



2P(s) 



eKb*{s). 



(Here we ignore the initial condition contribution which 
is justified by the stability of the system.) In terms of 
quadratures 6^ = ^{b + b*) and b^ = ^{b- b*), we find 
that 

&Sut(^) = SSpa(^)&^(^), by^.is) = S^pa(^)&^(s) 
where (in agreement with @, eq. (7.2.26)]) 



B. Cavity 

In a rotating reference frame, a model for a de- 
tuned cavity is characterized by the parameters Gcav = 
(1, ^a,uja*a); i.e. fl_ = w, ri+ = 0, C_ = C+ = 0, 
5" = /. 

This corresponds to an annihilation-form system 



bout = ^/ja + bu 



(44) 



when driven by vacuum input b. The transfer function 
for this system may readily be computed to be 



■^cav,-(s) — I 7 I ■ ' 
S + ^ + lU! 



(45) 



which in doublcd-up form is 



'^cav('^) 



S~\--^-\-iLJ 









(46) 



Thus this system is 

Gcav = (/, A(V7/,0),-iA(iu;,0)) € £"^(1). 

C. Degenerate Parametric Amplifier 

We consider the model for a degenerate parametric am- 
plifier (DPA) sec. 7.2], which corresponds to a single 
oscillator G coupled to a single field with S = I, lo^ — 0, 



"Dpa(*) 



"dpa('*) 



The DPA can be implemented in a single-ended cavity 
and a case that is of our main interest in this paper is 
the idealized one (for a full discussion sec 0, 10.2.1.g]) 
where k, e — > cx) (in practice to be taken large) such that 
the ratio is constant. Rescaling k = kn^ and e = fceo is 
equivalent to replacing k by kq and e by eg and rescaling 



^dpa(s, k = fc^o, e = fceo) 



2dpa(t, ko, Co)- 
k 



The limit fc — > cx) is the appropriate limit and here the 
cavity has an instantaneous response and the internal 
cavity dynamics are essentially eliminated by adiabatic 
elimination. This results in 6out(s) being given as the 
following Bogoliubov transformation of the input: 



bh{s) — — cosh(ro) b{s) — sinh(ro) b\s), 



where 



ro = In 



^0 + ep 



The output is then an ideal squeezed white noise process 
satisfying the quantum Ito rule discussed in section llll A| 
where 



N = sinh^ ro 
M = cosh To sinh rg = 



4^060 

2^060(^^0 — Cq) 
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Note here that M and A'' satisfy the relation |Mp = 
N{N + 1). In this hmit the DPA device behaves hke 
a static device that instantaneously outputs a squeezed 
white noise field from a vacuum white noise field source, 
and the transfer function has a constant Bogoliubov ma- 
trix value across all frequencies. That is 

~ - s 

5DPAstatic(s) = lim Sdpa(ti ^0, Co) 

fc^oo K 

= — A(coshro, sinhro), Vs G C, 
and the quadrature transfer functions are 

■^DPA static (•s) = ~e °, ■^DPA static (•^) ^ " ■ 

For a DPA device with a sufficiently wide bandwidth, 
one may approximately model it as a static DPA device 
with the ideal characteristics described above. Clearly, 
(n) is not closed with respect to this type of approx- 
imation. 



V. COMPONENTS INVOLVING BOGOLIUBOV 
TRANSFORMATIONS 

In section ITV CI we obtained the constant transfer func- 
tion Sdpa static G Sp(C) for a static approximation to a 
DPA. Such static approximations afford useful simplifica- 
tions, though in reality the DPA is a dynamical physical 
device. The idealized DPA therefore yields outputs that 
are a squeezing of the inputs. 

Motivated by this, in section IV Al we consider Bogoli- 
ubov matrices acting on boson fields, thereby extending 
the class of static components beyond unitary scattering 
devices. These components will be combined with linear 
dynamics in section IV CI to form a general class of quan- 
tum linear systems; such models may be useful when the 
time scales of the dynamical parts are slower than the 
time scales of the systems represented by static Bogoli- 
ubov matrices. These components will be combined with 
linear dynamics in section IV CI to form a general class 
of quantum linear systems; such models may be useful 
when the time scales of the dynamical parts are slower 
than the time scales of the systems represented by static 
Bogoliubov matrices. 



A. Bogoliubov Static Components 

More generally we could consider a static component 
which performs a Bogoliubov transformation of the input 
field <E !F{n): 

knt{t) ^ Sh^{t). (47) 

where now S G Sp(C"). This transformation, of course, 
preserves the canonical commutation relations so that 
bout G ^(n). 



Some caution should be applied here as we are now 
using the symbol S in (j47p in a purely algebraic man- 
ner as an element of Sp(C") when we strictly mean the 
second-quantization of the Bogoliubov matrix as an op- 
erator on the fields. Despite its formal similarity to ([6]), 
the relation (|T7)) is of a different character as the fields 
carry a continuous time-variable. Moreover, since such 
a transformation in general form linear combinations of 
field annihilation operator and creation operators, the 
transformation &out(0 = 5'&in(i) cannot be described by 
the usual Hudson-Parthasarathy QSDE for open Markov 
systems (cf. Section lTlI Cp . Such a QSDE can only model 
linear combinations of field annihilation operators of the 
form boutit) = A(5, 0)6in(i) for a unitary matrix S that 
appears as one of the parameters of the QSDE (here we 
set the other parameters to L = and H = Q). As such, 
in the transformation of fields with a non- unitary Bogoli- 
ubov matrix we do not have an analogue of ^ in the form 
of bout{t) = U {t)*bin{t)U{t) for some unitary process U{t) 
on the system and noise Hilbert space. At present we do 
not know whether a unitary transformation exists and if 
it exists what kind of dynamical equations it would sat- 
isfy. Since unitary evolution is a fundamental postulate 
of quantum mechanics, the situation is somewhat unsat- 
isfactory and is the subject of continuing research. How- 
ever, the relation (|17|) is nevertheless a useful idealization 
for certain devices used in quantum optics, such as what 
we have seen with the static DPA in Section IIV CI and 
has formally been employed upto now (see, for instance, 
the discussion in Chapter 7 of [6| on various quantum op- 
tical amplifiers). The physical meaning of the Bogoliubov 
transformation (j47p is correctly interpreted as a limiting 
situation. 



B. Bogoliubov Static Components as Limits of 
Dynamical Components 

The class of linear dynamical components described 
in section fill CI is not closed under input-output conver- 
gence. We now show how arbitrary static Bogoliubov 
components may arise as limits of unitary models. The 
idea is to exploit the Shale decomposition ([7]). Thus 
any given Bogoliubov matrix S has the decomposition 
S = A(5't^t, 0)A(coshi?, sinhi?)A(5in, 0), where Sin and 
S'out are some unitary matrices and R is some real diag- 
onal matrix. We note that the end terms A(S'j[jjj, 0) and 
A(S'in,0) can each be realized as a static passive net- 
work made of beamsplitters, mirrors and phase shifters. 
The middle term of course describes squeezing but this 
arises from a straightforward construction involving n in- 
dependent static-limit DPAs acting as ideal squeezing de- 
vices. (Each DPA corresponding to a diagonal entry of 
R providing a degree of squeezing (cf. Section IIV C\i as 
determined by that entry.) Then we note that we may 
approximate each DPA with a corresponding dynamic 
(non-ideal) DPA with appropriate parameters (sec the 
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discussion of the DPA in Subsection IIV Cp . 



C. Dynamical Bogoliubov Components 

We introduce an extension of the class of dynamical 
linear models C^^{n) considered up to now to accom- 
modate the notion of squeezing. This extension is in- 
spired by the factorization (j35p for open linear systems 
of (n) type, suggesting that we consider a new class 
of dynamical components of the form 



(48) 



where iI,C,n) G C^^ {n,m) and S = A{S-,S+) G 
Sp(C")- 

A system G ~ {S,C,Cl) is defined by equations (pS)) . 
where A, B and C are as before (section [ill Cp . but now 
D = S. We use the notation £^°s {n,m) to denote this 
class of systems and write C^°s {n) = UmC^°^' {n,m). 
The class C^°s {n) includes C^^ (n) special case 

(with S — A{S,0)). The justification for the cascade 
expression (j48p will be given in section IVIFi where we 
shall extend the series product for cascaded systems in 
£B°g- (n). 

The doubled up input-output map is of the form (|36p . 
where now D = S. The transfer function is explicitly 



Sg(s) 



A 




C 


s 



{s)^^C{sI-A)-^C^S + S. 

(49) 



The transfer function has the following properties: 



1. S 



G 



'A 






I 



2. Whenever its value exists, we have Sg(iw) € 
Sp(C"), for cj € K. 



Property 1 follows directly from ([49]) while Prop- 
erty 2 follows mutatis mutandis from the proof of |15l . 
Lemma 2] by the replacing with ^ and (A, S, C, D) with 
(i,B,(7,D). 

Physically, the meaning of ([48]) is that to process an 
input signal, the Bogoliubov transformation S is applied 
and the result is fed into the dynamical subsystem. As 
we have argued in section [VAl Shale's theorem precludes 
a unitary stochastic dynamical model giving rise to a 
system G G C^°k- {n). Nevertheless, as we have also seen 
that there will exist a sequences Gl e (n) such that 
pointwise 

with G e £^°s- (n) but not in C^^ (n). One might en- 
visage other modes of convergence of transfer function, 
however, we shall restrict to pointwise convergence for 
the purposes of this paper. It is interesting to note that 



the class of Hudson-Parthasarathy models is not closed in 
the above sense of convergence in the input-output sense, 
but may be extended to include Bogoliubov transforma- 
tions. 

We remark that in many cases where a boson field is 
in a squeezed state (recall section [III Bp . this field may 
be regarded as the output of a static Bogoliubov com- 
ponent S driven by vacuum inputs. This means, for ex- 
ample, that a dynamical component (1,C,J7) driven by 
squeezed fields may be represented as a dynamical Bo- 
goliubov component {S, C, fl). 



D. Example: Cavity with Squeezed Input 



Consider the cavity discussed in section lIVBi where 
now we suppose that the cavity input is given by the 
output of a squeezer Ggq, described by the Bogoliubov 
transformation 



A(cosh A, sinh A) 



cosh r sinh r 
sinhr coshr 



(50) 



That is, 



Gsq = (A(cosh A, sinh A), 0, 0) . 

The squeezed-input cavity Gcav, sq = Gc 
transfer function 



< Gsq has 



"^cav ('5)*5'sq 
s—^-\-iuj 

S — -^—lUJ 



coshr 
sinhr 



sinhr 
cosh r 



.(51) 



This corresponds to the equations 



a 




a* 












-(i 




iuj) 


cosh r sinh r 




■ h ' 


sinh r cosh r 




b* 





a 




a* 



bout 











a 


h* 
. out 







V7. 




a* 



+ 



cosh r sinh r 
sinh r cosh r 





■ b ' 




b* 



(52) 



The physical parameters for the squeezed-input cavity 

: Gcav,sq+ — 



are ^cav,sq- ^1 ^cav,sq+ 0; Gcav,sq 

0, *S'rav.S( 



= Ssq ~ A(coshr, sinhr), and so 
(A(cosh A, sinh A), 0), -iA{ioj, 0)). 

This system is a member of £^°s (l) but not of £"^(1). 



G, 



cav, sq 



VI. LINEAR QUANTUM FEEDBACK 
NETWORKS 



We are now in a position to described feedback net- 
works constructed from Bogoliubov dynamical compo- 
nents as nodes, and boson fields as links. The general 
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form of such a linear quantum feedback network (LQFN) 
is shown in Figure [D The fundamental algebraic tool for 
describing such networks in subsequent sections is the 
linear fractional transformation (LFT). 




B. Fractional Linear Transformations 

In this section we provide some technical results needed 
for the network theory described in subsequent sections. 

Lemma 1 Let S = A{S-,S+) G Sp(C"i+"^) with block 
decomposition 



where SJ^^ e C".x"fc. 
C2n,x2nfc^ have that 



'Setting Sjk = A{S^^,SJ'^) G 



SkiSkj = ^ SikS 



(53) 



k=l,2 



k=l,2 



FIG. 2; General form of a linear quantum feedback network 
(LQFN), with the time delay due to the spatial extent of the 
feedback connection represented by (see text). 



A. Linear Fractional Transformations 

Linear fractional transformations (LFTs) arise natu- 
rally when dealing with feedback networks, and a formal 
notation has been developed in classical linear systems 
theory, [23|. Consider a classical transfer function S(s) 
partitioned as 



S(.) 



Sii(s) Si2(.s) 
S2l(s) S22(s) 



corresponding to a partition u = (^1,^2)^, y = (2/1,2/2)^ 
of the input and output signals. If the system is placed in 
a feedback arrangement defined by U2 = K{s)y2, then the 
closed-loop system is described by the transfer function 

Enis) + E,2{s)K{s)[L ~ E22{s)Kis)]-^E2i{s). 

The arrangement is said to be well-posed whenever the 
inverse [L — 'E,22{s)K{s)]~^ exist. This transfer function 
is obtained by eliminating the in-loop variables. 

In what follows wc generalize this type of representa- 
tion to our class of LQFNs (see also Q, [111)- 



Proof. The relation S^S = L implies that 
A{Sl,-Sl)A{S-,S+) = A(/,0), and so SlS- - 
SjS* = I, SiS+ - SjS* = 0. These may be written as 

Y^k=l,2i^ki' ^kj ^ ^kJ ^kf") ~ J2k=l, 2i^ ki ^kj ^ 

S+^S^f) = 0. Therefore 



k=l,2 



k=l,2 



E ^(^ki^kj ^ki^ ^kf" ^ ^ki ^kj ^kJ^kf") 



fc=l,2 



which equals A ((5^ ,0) = Sij. The second identity simi- 
larly follows from SS^ = /. □ 

Theorem 2 Let S e Sp(C"i+"2) and define the 
fractional linear (Mobius) transformation : 
dom(«'|-*i) G C"2^"2 1-^ i)y 



- + Si2X{L - S22Xy'S2i, (54) 

with X € dom{'^^£^^) if and only if the inverse (1 — 
S22Xy^ exists. r/ien*|-"i maps Sp{C''^)ndom{^'f^^) 
into Sp(C"i). 

Proof. Wc first note the Siegcl-typc identities 



^f^\xf^^f\Y) = L-Sl^{L-X^Sl2)~\l-X^Y){L~S22Y)-^S2i, 
^Y^W^T^i^f = L - Si2{I - XS22y^{I - XY^){I - Sl2Y^y^S\2- (55) 

I 



These are structurally the same as the standard Siegel volution b replacing the usual Hermitian involution f , 
identities based on partitioning a unitary S but the in- 
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see Theorem 21.16 and Corollary 21.17 of reference 
The identities rely on the unitary analogue of the iden- 
tities (|53|) and so follow mutatis mutandis. Evidently, if 
X e Sp(C"=) ndom(*|-"i) then ^~'\Xf^-^\X) = 

^-^\x)^-^^xf = I. □ 



Corollary 3 IfK{iuj) is an Sp{C"^) -valued transfer ma- 
trix function taking values in c?om(^'|^^) for all to real, 

then the fractional linear transformation ^"^^^ [K (iuj)) 



will be Sp(C"i 



valued function of to. 

2-'!^ then 



In particular, if I (z dom{'^ - 



= Su + Su{l - S22)-'S2i e Sp(C"0- (56) 



C. Finite Time-Delay LQFNs 

A general LQFN is a network of linear quantum com- 
ponents Gv e £^°^ {n), labeled by the vertices v of the 
network, with quantum fields traveling along the edges. 
The edges are directed so that we distinguish inputs and 
outputs, and the multiplicity of input fields equals the 
multiplicity of outputs for each component. 

In a physical LQFN we will have time delays associ- 
ated with each internal edge due to the finite time taken 
by light to travel from an output port to an input port. 
In fact, we may lump the individual components as one 
single global component G with all external inputs going 
into a collective input port 1 and coming out from a col- 
lective output port 1, as in Figure [H Likewise, all the 
internal fields can be viewed as traveling from the col- 
lective output port 2 to the collective input port 2. The 
effect of the (multichannel) time-delay can be described 
by the operator O,- defined by 

e.(/l(t), . . . , fn{t)V = {flit - n), . . . , fn{t - r„))^, 

where ti > 0, . . . , t„ > are the time delays of each 
channel. Here fk (t) denotes the quantum stochastic pro- 
cess propagating along channel k in doubled-up form. 
For instance, fk(t) could be yfc(t), the doubled- up output 
quantum output processes propagating along channel k. 
In a slight abuse of notation, wc also occasionally over- 
load the notation 8r to denote the delayed version of a 
quantum process that is not in doubled-up form, such 
as when fk{t) is taken to be yk{t) for all k. Note that 



[Qri'iuj)]jk = 



Ojk- 



Extending the standard notation 



recalled in section fVI Ap . we denote this by 

Nr^d{G,er)- 

A Hamiltonian for a LQFN with squeezing compo- 
nents could be constructed approximately by replacing 
Bogoliubov components S with dynamical components 
G|. This would then fit into the QFN framework of d]. 



D. Parameters for Network Model 

We now suppose that the LQFN of Figure [2] is de- 
scribed by field channels 6i,6out,i and &2,&out,2 having 
lengths rii and n2, respectively, so that the total num- 
ber is 711+71,2 = n. The system is parameterized by 
G = {S, L, H), with S = A(S'_, 5"+) and we partition the 
matrices as 



C^ 



The field-field component of the input-output relations 
can be now written as 







^11 


•5 12 




. ^21 


•^22 . 



&out,i = ^ Gij{s)hj, 

3 = 1.2 



with transfer matrix function 



A 



Gi 

C2 



'[ClC\]S 



s 



(^) 



Ci 

C2 



{sI-A)-'[C\,C\]S + S 



(57) 



where 



S.k^HSjk.S^,), C,^^{Cj,Gp. 



The network Nr is given by the linear fractional trans- 
formation 



Ar 


G^Sr 


Gr 


Sr 



(58) 



where 



Sr = ^Il+5i2e,(/-522e0"'521 (59) 
Cr = Ci+Sr2Qr{I - S22Qry^C2 (60) 
Ar = A- C^jSj2er{I~S22er)-'C2. (61) 

Due to the non-zero delay, the network model Nr is non- 
Markovian. 



E. Zero Delay Limit Models 

Of particular interest are the simpler models that arise 
in the zero-delay limit 9t ^ / (t — > 0). Assume that 
I — Sii is invcrtible. From above we have 



iVo=m/)= (62) 



' Ao 


—CqSo 


Co 


So 



where 



^0 

Co 
Ao 



Sll+Si2{I-S22)-'S21 (63) 

Cl+^12(/-522)-'C'2 (64) 

^- E C'^S,2{I~S22)-'C2. (65) 

3 = 1,2 
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We note that 



An = 



where 



j=1.2 



(66) 



(67) 



Here, Im^X means 57 (^ — X^). The matrix Sq defined 
by (|63p is a Bogofiubov matrix, as it corresponds to the 
matrix in (|56p . Therefore the zero-delay Umit Nq is a 
Markovian system belongs to C^°s- (tj) with parameters 

-^0 = (^Sii + Si2{I — S22) ^821, 
Cl + Si2{I — S22) ^C2, 

+ Imt V C^jS,2{I - S22)-'C2) . (68) 

Thus £^°^ (ri) is closed with respect to this zero-delay 
limit network construction. 

Other types of limits are also considered in applications 
(see, e.g. [l|, sec. 2.3]). Suppose that the system G"^ and 
the delay r*^ depend on a small parameter e > 0, defining 
a physical regime of operation. Then one may obtain 
a limit model lim^^Q ^{C^ , Qt'^) ■ An example of this is 
considered in section fVII Bl below. 

F. Series Product 

The series product G2 < Gi of two systems Gi ~ 
[Si, Gi, fii) and Gi = (6*2, G2, ^2) follows from the zero- 
delay limit (p^ . For the series product, we interchange 
the index 1 and 2 (this simply means interchanging the 
role of Gi and G2 in the LQFN) and then set 5*12 = 5*1 , 
S'21 = S2, S22 = and Sii = (note that without 
the interchange we would be computing Gi < G2 instead 
G2<Gi). 

Substituting into (|63l64l65p we find 



(69) 



The matrices for G2 <i Gi = (5'sciics, Gsorios, ^^scrios) are 
given by 

S'scrios — S2S1, 

C'serics = G2 + S'2Gi = A(Gsorios-i C'soricsH-)) 



' A~ClS2Ci 


- 


(g^52 + Cl) 


Si' 


C'2 + S'2Gi 


S2S1 _ 



f^serics — ill + il2 + Im[,G2S'2Gi 
— ^(^seriGS— : ^series+) : 



where 



Gsorics- — G2- + S'2-Gi_ -|- S'2+G^, 
Gsorics+ = G2+ + S2-C1+ + S2+Cf_. 



From 

CIS2C1 = A(Gt,-C2\)A(52-,^2+)A(Gi_,Gi+) 
= A(X_,X+) 

with 

= {C\_S2^—C2+Sfj^)Ci^ + {C\_S2+~C2+Sf_)Cf^, 
we see that 

^^series- = ^^1- + ^2- + Tri^- ^ ^-): 

Succinctly, the series product in (n) is given by 

G2 <J Gi = (52^1, G2 + ^2Gi, fii +(12+ Imb{G^^2Gi}); 

(70) 

cf. HMD- 

Clearly, £^°s- (n) is a group with respect to the series 
product, and the classes of components £^^(n), Sp(C"), 
and U (n) are subgroups. If G = {S, G, fJ) the inverse is 
given by 



with transfer function given by 



(71) 
(72) 



The series product is not limited simply to feedforward, 
and the above formula apply to the case where the two 
systems have one or more modes in common. The series 
product is therefore highly non-trivial [1, [^, [l^ . 



G. Series Product and Cascaded Transfer 
Functions 

In classical linear systems theory, the transfer function 
of a cascade of two separate systems is obtained by mul- 
tiplying the transfer functions (see, e.g., [20|). 

We need to emphasize here that two systems will 
be distinct if they consist of different oscillator modes. 
Specifically let there be nii modes in the first system 

ai 



0-2 



and 1712 in the second, and set a = 
then we wish to consider Gi 

(S2,C2,n2) with 

A([Gi_,0],[Gi+,0]), 
-A 

A([0,G2 



e S (mi + 7712), 



(5i,Gi,fii) and G2 



Gi 
C2 









' 0' 




( 










) 



[0,G2+]), 



-in2 = -A 





■ ■ 




■ ■ 




( 


012- 


1 


012+ 


) 
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with respect to the decomposition C™ = C™i © C™^ 
This is simply a statement that the dynamics of Gi does 
not depend on the internal variables of G2, and vice- 
versa. Putting this particular form into (|69p . we then 
obtain 



Sg2<]Gi(s) 



where C, 



Ai 

'C2S2C1 A2 



-C'2S2Si 



S2CI, C2 
A{Cj-,Cj^+) and Aj 



S2S1 



(73) 



1 /^b/^ 



«A(r2j._, rij._|_). Here we use the unitary transformation 



ai 

0.2 



to present the transfer function in a more 



convenient form. The algebra is then similar to Section 
IV A in [Hi. 

If the systems are not separate, then we do not expect 
such a factorization of the transfer function to hold. The 
series product [1], [§], [l^ is defined quite generally in 
terms of physical parameters (section IIIIEp which may, 
for example, depend on the same oscillator mode vari- 
ables). In the general case, the transfer function can be 
computed using the general formulas ((62|) -([65 |) . 

Let us remark that the scries product inverse 
given by (|71[) may be realized in terms of a physical sys- 
tem that is not separate from the original system. Phys- 
ically, if we pass input fields through a system with pa- 
rameters C?, then gives the parameters required to 
undo the effect by passing the output back through the 
same system for a second pass. 



H. Inverse Transfer Functions 

The input-output relation 6out = 26in + S,a (0) may be 
inverted to yield 

6i„ = S-iL,t- 2-1^0(0). 

In particular, wc can give the following useful description 
of . The linear equations (|28p in the time domain may 
be rearranged algebraically to give 

j^a = [A- BD-^C)a + BD-%^t, 

with D invcrtiblc. and in the transform domain, wc de- 
duce that 



A 


B' 


-1 


'A-BD-^C 


BD-^' 


C 


D 




-D-^C 





For the model with parameters G = {S, C, (l), wc find 



'A 




-1 




-& 


c 






-s^c 





(s) ^ = I 



We note that ^g-i (*) — (^s) ^ = 



(74) 

Eq {s*)^ (recall 



1. Example: separate cavity inverse 

As a concrete example, consider the single mode cavity 
G ~ Gcav considered in IIVBI The transfer function ((45|) 
and related functions are given by 



Sg,-(s) 

^G-,-(.^) 



s — 


1 

2 


-f iuj 


s + 


7 
2 


-f iuj 


s ~ 


7 
2 


— iuj 


s + 


7 
2 


~ iuj 


s + 


7 
2 


+ iuj 


s — 


7 
2 


+ ILO 



That is, G "'^ is obtained from G by keeping C_ = .y/7, 
C+ = and replacing = a; by —uj. In what follows we 
obtain a physical realization G of the transfer function 
(s) that is a system distinct from G, so that 



(75) 



The system G is obtained by setting = — lj, fi^ = 
and swapping G_ = and G+ = 

In this example we can take the two modes to be 

G 5 (2) and write G as (^Si ~ I,Ci, Clij where 
Gi = A([V7,0],[0,0]), 





—iuj 




■ 


■ 




( 













) 



and G as ^82 = ^, G2, —^2j where 

G2 = A([0,0],[0,V^), 
-in2 = -A 





—iuj 




■ 


0' 




( 













) 



That is, in terms of scattering matrices, coupling opera- 
tors and Hamiltonians (1311) wc have 



G = (/, ^ai, wa*ai), and G = (/, ^a^, -^^0202). 

The series product G < G is given by 

I,C ^Ci+C2,n^ni+n2+ Imt[G^Gi]) (recall 

(|70|) ) and this corresponds to a system with a nontrivial 
dynamics, since (by (|34[) 



S^C{sI + A^)-'^C^ 



G<G ^ {I,y/^{ai + a2),uj{alai - 0302) +Im[a2ai]). 
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Nevertheless, some calculation shows that 

as required by (j73|) and (j75p . This continues to hold if 
there is a term f2_|_ = lu^ added to the cavity Hamilto- 
nian. 

2. General separate system inverses 



We now show how to obtain a physical realization for 
Sq^{s) in the case where = 0. We have 



= + s^c{si -^c^cy^c^ 



Now CC^ is not definite and in fact we have 

CC^ = A{c^cl~c+cl,-c^cl + c+cl) 

= -A(C+,C_)A(C|,-CT) 



where 



= A(C+,C_). 



(76) 



We therefore obtain 



1 



where 



Sp^s) = [I - S^{s + -KK^)-^KK^]S^ 

= [I- S^K{s + ^K^Ky^K^]S^ 
= 2g(s), 



G= {S\S^K,0). 



(77) 



We note that this choice of G is not unique in producting 
a transfer function inverse inverse to Sq- Finding an 
inverse when 57 7^ is more involved. 



3. Comments 

We have the mapping G 1-^ from the group of 
£^°s-(n) of systems parameters with series product to 
the group of matrix transfer functions. This mapping is 
not however a group homomorphism. Indeed, we typi- 
cally have 

Sg2<Gi(s) 7^ 2g2(s)Sgi(s), 

though equality - the cascade formula ([75)1 - holds when 
the systems are separate assemblies of oscillators. 

We should also caution that, as we have seen in the 
example in IVIH 11 there are solutions for G other than 
the trivial G ~ (/, 0, 0) to the equation (s) = /. 



VII. NETWORK EXAMPLES 

A. In-loop Squeezing and Cavity as a Feedback 
Network 

We now describe the LQFN of Figurefl] which contains 
a cavity and squeezer in a feedback loop resulting from 
interconnection with a beamsplitter. The total propaga- 
tion delay around the loop is r, which we take to be small 
and send t 0. In order to determine an equivalent 
zero-delay limit model, following the general approach of 
section rvTEl we re-draw the network as shown in Figure 

131 




FIG. 3: Linear quantum feedback network of Figure[T]redrawn 
in standard from as Figure (2] 



As indicated in Figure [31 the in-loop system G = 
Gdyn <i Ggtatic is 9, dynamical Bogoliubov component ob- 
tained by cascading the beamsplitter, the (augmented) 
squeezer (which together form Gstatic), and the (aug- 
mented) cavity Gdyn- 

The static part Gstatic is described as follows. Because 
the beamsplitter has two inputs and two outputs, we aug- 
ment the squeezer Sgq (given by (|50|) ) by including a di- 
rect feed through channel (vi to wi in Figure[3]). Because 
the squeezer is represented by a static Bogoliubov trans- 
formation expressed in doubled-up form, we express the 
beamsplitter in doubled-up form: St = A(S'fc,0). To be 
clear, the beamsplitter is described by 



Vl 
V2 



aui — (3u2 
Pui -f au2, 



where laP 



and 



1, a* (3 = P*a. Thus we have 



a -/3 
(3 a 



(78) 



(79) 



Sh — 



A(56,0) 

a — /3 

(5 a 








a* 






-P* 

a* 



(80) 



16 



The static component Gstatic has inputs {ui,U2)~^ and 
outputs {wi, ^2)^, and is given by the Bogohubov matrix 



R = A 

= A(i?_,i? 



1 

coshr 




sinh r 



(81) 



where 



R- = 



i?4 



a -(3 
(3 cosh r a cosh r 


(3* sinh r a* sinh r 



(82) 



The dynamic component Gdyn, with inputs (i«i,i(J2) 
and outputs (2/1, 2/2)^, is given by 



a 




■-(?+*^) 







a 


a* 











a* 



77 
^7 



Wl 

wl 

W*2 



Thus yl_ 



" yi ' 




" 


■ 




" Wl ' 






V7 







a 


+ 


W2 














a* 




wl 


-y*2- 











. ^2 . 


L_ = 




zo;), 




= 0, i = 


= A(-(i 



(83) 





x/7 



C4 



(5 = A( 





V7 



,0). (84) 



Also, Q,- = cj, fi+ = 0. 

Now that we have a complete model for the in-loop 
system G, we may apply the formulas in section fVI El to 
obtain a zero-delay network model Nq. This involves first 
working out the Bogoliubov matrix in partitioned form: 



= A(a,0), 5i2 = A(-/3,0), 
5*21 A(/3coshr,/3* sinhr), 
5*22 = A(a coshr, a* sinh r). 



(85) 



We set a — ^/e and (3 = — e to simplify some of the 
algebra. The network model is given as follows. We now 
use the equations (|63|64l65p to determine the network 
parameters. The equivalent network Bogoliubov matrix 
is 



,§0 — A(Ve— - — ^ (coshr 



Vi),-^i^^^), (86) 



where 

Next, 
Go 



1 — 2\/ecosh7 



-A(l — -^e cosh r, -^e sinh r). 



(87) 



so that 



Gn 



G+ 



vT 



(1 — Ve coshr). 



\/e sinh r. 



Now io = A(Ao , yl(j"), where 



Vf7 



From this we compute 



(cosh r — ^/e), Aq — sinh r 

(89) 



t smh r. 



(90) 



We therefore see that not only does the network model 
iVo G have a non-trivial static Bogohubov term, 

it also has field couplings involving a creation operator 
a*, and Hamiltonian terms involving and (a*)^. 

Stability of the feedback system may be analyzed using 
the methods of section IIII Dl or the small gain theorem 

As a possible application, we note that the squeezing 
parameter of a DPA may be altered by placing it in-loop 
in a beam-splitter arrangement of this type [25|. 



B. Dynamics from Feedback 

In this example we give an illustration from quantum 
optics showing that LQFNs involving only static com- 
ponents may give rise to dynamical behavior. This dy- 
namical behavior is due to a time delay in the feedback 
loop. We consider the network shown in Figure [1] p^ . 
0. This is a special case of the LQFN network of Figure 
[U but with no squeezing and no cavity. The beamsplitter 
56 is given by ^ or[7Hl with a = Ve and /3 = ^/l - e. 




FIG. 4: Beamsplitter feedback network with time delay r. 



Feedback in the network is defined by the constraint 
U2{t) = Qry'iit) ~ V2{t ^ t), where r > is the time 
taken for light to travel from the output to the input. 
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This network is a LQFN with = (A(S'6, 0), 0, 0) S 
U{2). With e > fixed, the zero-delay network model 
T ^ is the system Nq = di{G^ J) = (A(-/, 0), 0, 0) e 
U{1), with transfer function No{s) = — 1, a trivial pass- 
through system with sign change (phase shift). 

Now if the reflectivity coefficient e and the time delay 
arc comparable, say r = e/7, where 7 > 0, then we 
obtain a dynamical model as e — > 0, [l|. Section 2.3] (recall 
section IVIE] above). Indeed, solving (|79|) and ([80)) in the 
frequency domain we find that the transfer function is 



ee 



-sc/7 



1 - ^/l - ee-^'^/'' 
(91) 

By L'Hopital's rule, we find that the limit transfer func- 
tion is 

N{s) = IhnN^is) = 1 - —L^ = '—2/1. (92) 

This transfer function corresponds to a cavity N = 
(/,^a,0) = (1,^1,0) e /:"P(1,1), where a e 5(1), 
Q, [l|. Here, 7 plays the role of the coupling strength 
between the trapped cavity mode and the external free 
field. 

This example shows that U (2) is not closed under 
this type of physically natural approximation process 
(since the limit belongs to £^^(1) which is outside U{2)), 



while £^°s (2) is closed (since it contains both 17(2) and 

/:HP(i)). 



VIII. CONCLUSION 

We have shown how to extend linear quantum dynam- 
ical network theory to include static Bogoliubov compo- 
nents (such as squeezers). This unified framework ac- 
commodates squeezing components which are important 
in quantum information applications. We provided tools 
for describing network connections and feedback using 
generalizations of linear fractional transformations and 
the series product, 0, 0, [3, [10, Chapter 10]. We 
have also defined input output maps and transfer func- 
tions within this linear quantum network theory, and 
shown how they can be used in applications. Finally, 
we explained the natural group structure arising from 
the series product. 
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To see this, let S > be a complex 
Ell ~ 0. Set It (t) = {txi,X2, ■ ■ • , Xk) 



k X k matrix with 
then we have < 



u{tyEu{t) = 2tReJ2jyiXiEijXj + J2j,k>i^jEjkXk. 
However, the only way to guarantee this inequality for all 

^^j^ o^ij^ij^j — 0. Replacing 



real t is to require that ReX]j>i xiE 
t with it shows that xiEijXj 
As this must be true for arbitrary Xj, 
Eij = E*i = for all j > 1. 



likewise vanishes, 
we conclude that 



